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The preparative method, characterization and magnetic susceptibility \ measurements versus tem- 
perature T of the heavy fermion transition metal oxide LiV2 04 are reported in detail. The intrinsic 
x(T) shows a nearly T-independent behavior below ~ 30 K with a shallow broad maximum at 
pa 16 K, whereas Curie- Weiss-like behavior is observed above ~ 50-100 K. Field-cooled and zero- 
field-cooled magnetization Af obs measurements in applied magnetic fields H — 10-100 G from 1.8 
to 50 K showed no evidence for spin-glass ordering. Crystalline electric field theory for an assumed 
cubic V point group symmetry is found insufficient to describe the observed temperature variation 
of the effective magnetic moment. The Kondo and Coqblin-Schrieffer models do not describe the 
magnitude and T dependence of % with realistic parameters. In the high T range, fits of x(T) 
by the predictions of high temperature series expansion calculations provide estimates of the V-V 
antiferromagnetic exchange coupling constant J/ks ~ 20 K, g- factor g ~ 2 and the T-independent 
susceptibility. Other possible models to describe the x(T) are discussed. The paramagnetic impuri- 
ties in the samples were characterized using isothermal M obB (H) measurements with < H < 5.5 T 



at 2 to 6K. These impurities are inferred to have spin Simp 
concentrations of 0.01 to 0.8%, depending on the sample. 

PACS numbers: 71.28.+d, 75.20.Hr, 61.66. Fn, 75.40.Cx 
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I. INTRODUCTION 

Especially since the discoveries of heavy fermion (HFIl 
and high temperature superconducting compounds ,E 
strongly correlated electron systems have drawn much 
attention both theoretically and experimentally. Ex- 
tensive investigations have been done on many cerium- 
and uranium-based HF compounds a The term "heavy 
fermion" refers to the large quasiparticle effective mass 
m* /m c ~ 100-1000 of these compounds inferred from the 
electronic specific heat coefficient j(T) = C C (T) /T at low 
temperature T, where m e is the free electron mass and 
C c is the electronic specific heat. Fermi liquid (FL) the- 
ory explains well the low-T properties of many HF com- 
pounds. Non-FL compoundso are currently under interu 
sive study in relation to quantum critical phenomenal 
The transition jnetal oxide compound L1V2O4 was re- 
cently reportedQ to be the first d-electron metal to show 
heavy FL behaviors characteristic of those of the heaviest 
mass /-electron systems. 

L1V2O4 has the face-centered-cubic (fee), normal- 
spinel structure with space group FdSm [Fig. [y£a)], first 
synthesized by Reuter and Jaskowsky in 1960.0 The V 
ions have a formal oxidation state of +3.5, assuming 
that those of Li and O are +1 and —2, respectively, 
corresponding to 1.5 d-electrons per V ion. In the nor- 
mal oxide spinel L1V2O4, the oxygen ions constitute a 
nearly cpjbic-close-packed array. Lithium occupies the 
8a sitesja corresponding to one-eighth of the 64 tetrahe- 
dral holes formed by the close-packed oxygen sublattice 
in a Bravais unit cell that contains eight Li[v~2]04 for- 
mula units. Vanadium occupies the 16d sites (enclosed 
in square brackets in the formula) , corresponding to one- 
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FIG. 1. (color) (a) Normal spinel structure of L1V2O4 with 
an fee Bravais unit cell, (b) A part of the structure depict- 
ing the trigonally distorted oxygen octahedra. The distortion 
shown is exaggerated for clarity and corresponds to an oxy- 
gen parameter u = 0.27. Small, medium and large spheres 
represent lithium, vanadium and oxygen, respectively; their 
sizes have no intended physical significance. 
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half of the 32 octahedral holes in the oxygen sublattice 
per unit cell. All of the V ions are crystallographically 
equivalent. Due to this fact and the non-integral V ox- 
idation state, the compound is expected to be metal- 
lic, which was confirmed by single-crystal resistivity p(T) 
measurements by Rogers et alu The V atoms constitute 
a three-dimensional network of corner-shared tetrahedra. 
The L1V2 sublattice is identical to the cubic Laves phase 
(C15) structure, and the V sublattice is identical with 
the transition metal T sublattice of the fee -R2T2O7 py- 
rochlore structure. 

Despite its metallic character, L1V2O4 exhibits a 
strongly temperature dependent magnetic susceptibility, 
indicating strong electron correlations. In the work re- 
ported before 1997, the observed magnetic susceptibil- 
ity x° bs (^) was found to increase monotonically with 
decreasing T down to « .i^K-and to approximately-, fol- 
low the Curie- Weiss law.EJ~t3 Kessler and SienkoE2l in- 
terpreted their x° bs (T) data as the sum of a Curie- Weiss 
term 2C/(T — 0) and a temperature- independent term 
Xo = 0.4 x 10~ 4 cm 3 /mol. Their Curie constant C was 
0.468 cm 3 K/(molV), corresponding to a V +4 g-factor of 
2.23 with spin S = 1/2. The negative Weiss tempera- 
ture 6 — — 63 K suggests antiferromagnetic (AF) inter- 
actions between the V spins. However, no magnetic or- 
dering was found above 4.2 K. This may be understood 
in terms of possible suppression of long-range magnetic 
ordering due to the geometric frustration among,the AF- 
coupled V spins in the tetrahedra network. Ej'tZI Similar 
values of-Cpand have also been obtained by subsequent 
workersJiirEij as shown in Table |, in which reported crys- 
tallographic dataEjt 2 ] are also shown. This local mag- 
netic moment behavior of LiV2C>4 is in marked contrast 
to the magnetic properties of isostructural LiTi204 



which manifests a comparatively temperature indepen- 
dent Papk' paramagnetism and superconductivity (T c < 
13.7 K)B 

Strong electronj-cpcrelations in LiV204 were inferred 
by Fujimori et aZ.cjta from their ultraviolet (UPS) and 
x-ray (XPS) photoemission spectroscopy measurements. 
An anomalously small density of states at the Fermi 
level was observed at room temperature which they at- 
tributed to the effect of long-range Coulomb interactions. 
They interpreted the observed spectra assuming charge 
fluctuations between d 1 (V 4+ ) and d 2 (V 3+ ) configura- 
tions on a time scale longer than that of photoemission 
(~ 10~ 15 sec). Moreover, the intra-atomic Coulomb re- 
pulsion energy, U, was found to be ~ 2 eV. This value is 
close to the width W of the t2 g conduction band 

calculated for LiTi2 04.l 2 3eZl From these observations, one 
might infer that U ~ W for L1V2O4, suggesting proxim- 
ity to a metal-insulator transition. 

We and collaborators recently reported that LiV204 
samples with high magnetic purity display a crossover 
from the aforementioned localized moment behavior 
above ~ 100 K to a nearly, temperature independent sus- 
ceptibility below ~ 3OK.0 This new finding was also re- 
ported independently and nearly simultaneously by two 
other groups. EHlEa Specific heat measurements revealed a 
rapidly increasing j(T) with decreasing temperature be- 
low ~ 30 K with an exceptionally large value 7(1 K) »s 
0.42 J/molK 2 .l3 To our knowledge, this 7(1 K) is the 
largest value reported for any metallic d-electron com- 
pound, e.g., Y .9 7 Sc .o 3 Mn2 (<J12J/mol K 2 )(Refj2|) 
and V 2 -y0 3 (< 0.07 J/mol K 2 )EJ The Wilson ratioSJat 
low T was found to be i?w ~ 1.7, consistent with a heavy 
FL interpretation. From 7 Li NMR measurements, the T 
variation of the Knight shift K was found to approxi- 



TABLE I. Lattice parameter 00, oxygen parameter u (see text) and magnetic parameters xo, C and 9 reported in the 
literature for LiV2C>4. The At values shown are for the second setting of the space group Fd3m from the International Tables 
for Crystallography, Vol. A.U The "T range" is the temperature range over which the fits to the susceptibility data were done, 
Xo is the temperature-independent contribution, C is the Curie constant and 9 is the Weiss temperature. The error in the last 
digit of a quantity is given in parentheses. Unless otherwise noted, all measurements were done on poly crystalline samples. 
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a This value was digitized from the published figure. 

b Single crystal susceptibility data, corrected for the contribution of 10% V4O7. 
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mately follow that of the susceptibilityflli&il The 7 Li 
nuclear spin-lattice relaxation rate 1/Ti in LiV204 was 
found to be proportional to T below ~ 4K, with a Ko- 
rringa ra£i0-.Qa. the order of unity, again indicating FL 
behavior .&EJ~E3 

In this paper we present a detailed study of the syn- 
thesis, characterization and magnetic susceptibility of 
L1V2O4. In Sec. [njour synthesis method and other exper- 
imental techniques are described. Experi mental results 
and analyses are given in Sec. [II. In Sec. Ill A , after a 
brief overview of the spinel structure, we present struc- 
tural characterizations of nine L1V2O4 samples that were 
prepared in slightly different ways, based upon our results 
of thermogravimetric analysis (TGA), x-ray diffr action 
measurements and their Rictvcld analyses. In Sec. IIIB , 
results and anal yses of magnetization measurements are 
given. In Sec. [II B 1 an overview of the x° bs (T) = 
M° hs (T) I H d ata of a ll nine samples studied is presented. 
Then, in Sec. fn B 2\ we determine the magnetic impu- 
rity concentrations from analysis of the M° hs (H) data. 
Low-field (H = 10-100 G) x° hs {T) susceptibility data, 
measured after zero-field c ooling ( ZFC) and field cooling 
(FC), are presented in Sec. IIIB 3 ^, from which we infer 
that spin-glass ordering does not occur above 2K. The 
above determinations of magnetic impurity contributions 
to M° hs (H,T) allow us to extract the intrinsic suscep ti- 
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bility xCO from x° bs (T), as explained in Sec. 
The paramagnetic orbital Van Vl eck su sceptibility x x 
contribution is determined in Sec. IV A from-a^o-called 



K-x analysis using 51 V NMR measurements.EjS We at 
tempt to interpret the x(T) data using three theories. 
First, the predictions of high temperature series expan- 
sion (HTSE) calculations for the spin 5=1/2 Heisen- 
berg model are compared to our x(T) data in Sec. IV B . 
Second, a crystalline electric field theory prediction with 
the assumption of c ubic p oint symmetry of the vanadium 
ion is tested in Sec. IVC Third, we test the applicability 
of the Kond o and Coqblin-Schrieffer models to our x(T) 
data in Sec. IV D . A summary and discussion are given 
in Sec. [y| Throughout this paper, a "mol" means a mole 



of L1V2O4 formula units, unless otherwise noted. 



II. SYNTHESIS AND EXPERIMENTAL DETAILS 

Polycrystalline samples of L1V2O4 were prepared us- 
ing conventional solid-state reaction techniques with two 
slightly different paths to the products. The five samples 
used in our previous workD (samples 1 through 5) were 
prepared by the method in Ref. |23|. Two additional sam- 
ples (sampleSi-61 and 7) were synthesized by the method 
of Ueda et alB2 Different precursors are used in the two 
methods: "Li 2 V03.5" (see below) and Li 3 VC>4, respec- 
tively. Both methods successfully yielded high quality 
LiV2C>4 samples which showed the broad peak in x obs (T) 
at k 16 K. In this report, only the first synthesis method 
is explained in detail, and the reader is referred to Ref. |2^ 
for details of the second method. 

The starting materials were Li 2 C03 (99.999%, John- 
son Matthey), V 2 3 , and V 2 5 (99.995%, Johnson 
Matthcy). Oxygen vacancies tend to be present in com- 
mercially obtained V2O5E3 Therefore, the V2O5 was 
heated in an oxygen stream at 500-550 °C in order to fully 
oxidize and also dry it. V2O3 was made by reduction of 
either V 2 5 or NH4VO3 (99.995%, Johnson Matthey) in 
a tube furnace under 5 % H2 /95 % He gas flow. The heat- 
ing was done in two steps: at 635 °C for ss 1 day and then 
at 900-1000 °C for up to 3 days. The oxygen content of 
the nominal V2-yC>3 obtained was then determined by 
thermogravimeteric analysis (TGA, see below). The pre- 
cursor "Li2V0 3 .5" (found to be a mixture of LisV04 and 
LiV03 from an x-ray diffraction measurement) was pre- 
pared by heating a mixture of Li2C0 3 and V2O5 in a 
tube furnace under an oxygen stream at w 525 °C until 
the expected weight decrease due to the loss of carbon 
dioxide was obtained. Ideally the molar ratio of Li2C0 3 
to V2O5 for the nominal composition Li2V03.5 is 2 to 
1. A slight adjustment was, however, made to this ratio 
according to the actual measured oxygen content of the 
V 2 - y 3 (y ~ 0.005 to 0.017) so that the final product 



TABLE II. Results of Rietveld refinements of x-ray diffraction measurements and magnetization M ohs (H) isotherm analyses. 
The oxKgen parameter (u) is for the second setting of the space group Fd3m from the International Tables for Crystallography, 
Vol. A.a / str imp is the impurity concentration. The error in the last digit of a quantity is given in parentheses. The detection 
limit of /str imp is assumed to be 1%.l3 For samples 3 and 7 in which no discernable impurities were seen, this detection limit 
is listed; the Rietveld refinement for sample 5 directly yielded / str i mp < 1%. 
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is stoichiometric LiV204. This precursor and V2-j,03 
were ground thoroughly inside a helium-filled glovebox. 
The mixture was then pelletized, wrapped in a piece of 
gold foil, sealed into a quartz tube under vacuum, and 
heated between 570 °C and 700 °C for < 2 weeks. The as- 
prepared samples were all removed from the oven at the 
final furnace temperature and air-cooled to room tem- 
perature. For samples 2 and 3 additional heating at a 
higher T = 750 °C was given, with a repeated sequence of 
grinding, repelletizing and reheating for sample 2. From 
« 725 °C different methods of cooling, liquid-nitrogen or 
ice-water quenching or slow-oven cooling, were applied 
to pieces from sample 2, yielding samples 4, 4 A and 4B, 
respectively. 

Using a Rigaku Geigerflex diffractometer with a curved 
graphite crystal monochrometer, x-ray diffraction pat- 
terns were obtained at room temperature with Cu Ka 
radiation. Rietveld analyses of the diffraction patterns 
were carried out using the angle-dispersive jjj-rpy diffrac- 
tion version of the RIETAN-97/3 programOH 

TGA measurements were done using a Perkin-Elmer 
TGA 7 Thermogravimetric Analyzer. Oxygen contents of 
the samples were calculated from weight gains after heat- 
ing in an oxygen flow to 540 °C for L1V2O4 and 620 °C for 
Y2-yOs, assuming that the oxidized products contained 
vanadium as V +5 . 

Magnetization M obs measurements were performed us- 
ing a Quantum Design MPMS5 superconducting quan- 
tum interference device (SQUID) magnetometer over the 
T range from 1.8-2 K to 400 K with H up to 5.5 T. Zero- 
field-cooled (ZFC, usually obtained by quenching the su- 
perconducting solenoid) M ohs (H = 1T,T) scans were 
carried out and isothermal M° hs (H ) data at various tem- 
peratures were obtained. Low-field (10-100 G) ZFC and 
field-cooled (FC) M(T) scans were done from 1.8-2 K 
to 50 K in order to check for the presence or absence of 
spin-glass ordering. 



III. RESULTS AND ANALYSES 
A. Structure 

X-ray diffraction patterns of our nine L1V2O4 sam- 
ples revealed that the samples were single-phase or very 
nearly so. Figure ^(a) shows the diffraction pattern of 
sample 7 which has no detectable impurities. The nine 
samples described in detail in this paper are categorized 
into three groups in terms of purity: essentially impurity- 
free (samples 3 and 7), V3O5 impurity (samples 1, 4 and 
6) and V2O3 impurity (samples 2, 4A, 4B and 5). The 
presence of these impurity phases is detected in magnified 
views of the diffraction patterns as shown in Fig. ||(b). 
Results from Rietveld analyses of the diffraction patterns 
for these samples are given in Table |l| The refinements 
of the spinel phase (space group Fd3m, No. 227) were 



based on the assumption of exact L1V2O4 stoichiomc- 
try and the normal-spinel structure cation distribution. 
The values of the isotropic thermal-displacement param- 
eters B of lithium and oxygen were taken from the Ri- 
etveld analysis of neutron diffraction measurements on 
our L1V2O4 sample 5 by Chmaissem et al.,t3 and fixed 
throughout to £?Li = LI A and Bo — 0.48 A, respec- 
tively. These two atoms do not scatter x-rays strongly 
enough to allow accurate determinations of the B values 
from Rietveld refinements of our x-ray diffraction data. 
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FIG. 2. (a) X-ray diffraction pattern of L1V2O4 sample 7. 
The spinel-phase peaks are indexed as shown, (b) Expanded 
plots of the X-ray patterns of samples 1 (top), 2 (middle) 
and 7 (bottom). Indexed peaks are those of the spinel phase. 
Sample 1 has V3O5 impurity (filled circles), whereas sample 2 
has V2O3 impurity (filled squares). Sample 7 has no impurity 
peaks except possibly the very weak unidentified one marked 
with a star. 

The positions of the oxygen atoms within the unit cell 
of the spinel structure are described by a variable oxygen 
parameter u associated with the 32e positions in space 
group FdZm. The value of u [in the space group setting 
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with the origin at center (3m)] for each of our samples 
was found to be larger than the ideal close-packed-oxygen 
value of 1/4. Compared to the "ideal" structure with 
u = 1/4, the volumes of an oxygen tetrahedron and an 
octahedron become larger and smaller, respectively. The 
increase of the tetrahedron volume takes place in such 
a way that each of the four Li-0 bonds are lengthened 
along one of the <111> directions, so that the tetrahe- 
dron remains undistorted. As a result of this elongation, 
the tetrahedral and-pctahedral holes become respectively 
larger and smaller .c3 Each of the oxygen atoms in a tetra- 
hedron is also bonded to three V atoms. Since the frac- 
tional coordinates of both Li and V are fixed in terms of 
the unit cell edge, an oxygen octahedron centered by a V 
atom is accordingly trigonally distorted. This distortion 
is illustrated in Fig. 0(b). 

The nine LiV204 samples were given three different 
heat treatments after heating to 700 to 750 °C: air- 
cooling (samples 1, 2, 3, 5, 6 and 7), liquid-nitrogen 
quenching (sample 4), ice- water quenching (sample 4A) 
or oven-slow cooling at w 20 °C/hr (sample 4B). Possible 
loss of Li at the high synthesis temperature, perhaps in 
the form of a lithium oxide, was a concern—In a detailed 
neutron diffraction study, Dalton et aLcil determined 
the lithium contents in their samples of Lii +:E Ti2- 2 ;04 
(0 < x < 0.33), and found lithium deficiency in the 8a 
site of the spinel phase of all four samples studied. If 
the spinel phase in the Li-V-0 system is similarly Li- 
deficient, then samples of exact stoichiometry L1V2O4 
would contain V-0 impurity phase(s), which might then 
explain the presence of small amounts of V2O3 or V3O5 
impurity phases in most of our samples. 

Sample 3 was intentionally made slightly off- stoichio- 
metric, with the nominal composition LiV1.92O3.89. A 
TGA measurement in oxygen showed a weight gain of 
12.804% to the maximally oxidized state. If one as- 
sumes an actual initial composition LiVi. 9203.89+5, this 
weight gain corresponds to 8 — 0.08 and an actual ini- 
tial composition of LiV1.92O3.97 which can be rewritten 
as Li1.01V1.93O4 assuming no oxygen vacancies on the 
oxygen sublattice. On the other hand, if one assumes 
an actual initial composition of Li1_2V1.92O3.89, then 
the weight gain yields x — 0.19, and an initial composi- 
tion Lio.8iV1.92O3. 89 which can be similarly rewritten as 
Lio.83V1.97O4. Our Rietveld refinements could not dis- 
tinguish these possibilities from the stoichiometric com- 
position Li[V2]04 for the spinel phase. 

Sample 4, which was given a liquid-nitrogen quench 
from the final heating temperature of ~ 725 °C (la- 
belled "LN2" in Table Hj), is one of the structurally least 
pure samples (see Table |l]). Our Rietveld refinement of 
the x-ray diffraction pattern for this sample did not re- 
veal any discernable deviation of the cation occupancy 
from that of ideal Li[V2]04. There is a strong simi- 
larity among samples 4, 4A (ice-water quenched) and 
4B (oven-slow cooled), despite their different heat treat- 
ments. These samples all have much larger lattice pa- 
rameters (do > 8.246 A) than the other samples. The 



as-prepared sample 2, from which all three samples 4, 
4A and 4B were obtained by the above quenching heat 
treatments, has a much smaller lattice parameter. On 
the other hand, the oxygen parameters u of these four 
samples are similar to each other and to those of the 
other samples in Table ||. 

The weight gains on oxidizing our samples in oxygen 
in the TGA can be converted to values of the average 
oxidation state per vanadium atom, assuming the ideal 
stoichiometry L1V2O4 for the initial composition. The 
values, to an accuracy of ±0.01, are 3.57, 3.55, 3.60, 
3.56, 3.56, 3.57, 3.57, 3.55 for samples 1-7 and 4B, re- 
spectively. This measurement was not done for sample 
4A. These values are systematically higher than the ex- 
pected value of 3.50, possibly because the samples were 
not completely oxidized. Indeed, the oxidized products 
were gray-black, and upon crushing were brown, rather 
than a light color. On the other hand, x-ray diffraction 
patterns of the "L1V2O5.5" oxidation products showed 
only a mixture of L1VO3 and Li4Vio027 phases _as ex- 
pected from the known Li2 0-V2 0s phase diagramo Our 
upper temperature limit (540°C) during oxidation of the 
LiV204 samples was chosen to be low enough so that the 
oxidized product at that temperature contained no liq- 
uid phase; this temperature may have been too low for 
complete oxidation to occur. In contrast, our V2-j,03 
starting materials turned orange on oxidation, which is 
the same color as the V2O5 from which they were made 
by hydrogen reduction. 

B. Magnetization Measurements 

1. Overview of Observed Magnetic Susceptibility 

An overview of the observed ZFC magnetic suscepti- 
bilities x° hs (T) = M ohs (T)/H aX H = LOT from 1.8 
2K to 400 K of the nine L1V2O4 samples is shown in 
Figs. I (a), (b) and (c). The x° hs (T) data for the vari- 
ous samples show very similar Curie- Weiss-like behavior 
for T > 50 K. Differences in X (T) between the sam- 
ples appear at lower T, where variable Curie-like Ci mp /T 
upturns occur. 

Samples 1 and 6 clearly exhibit shallow broad peaks 
in x obs at T w 16 K. The x° hs (T) of sample 6 is system- 
atically slightly larger than that of sample 1; the reason 
for this shift is not known. Samples 3 and 4 also show 
the broad peak with a relatively small Curie-like upturn. 
Samples 2 and 7 show some evidence of the broad peak 
but the peak is partially masked by the upturn. For sam- 
ples 4 A, 4B and 5, the broad peak is evidently masked 
by larger Curie impurity contributions. From Fig. || and 
Table |l| the samples 1, 4 and 6 with the smallest Curie- 
like magnetic impurity contributions contain V3O5 impu- 
rities, whereas the other samples, with larger magnetic 
impurity contributions, contain V2O3 impurities. The 
reason for this correlation is not clear. The presence of 
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the vanadium oxide impurities by itself should not be a 
direct cause of the Curie-like upturns. The susceptibility 
of pure V2O3 follows the Curie- Weiss law in the metallic 
T region above ~ 170 K, but for T < 170 K it becomes 
an antifcrromagnetic insulator, showing a decrease in 
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FIG. 3. Observed magnetic susceptibility x° bB (T) 

(= M obs /H) of all the nine samples studied, measured with 
H = 1 T after being zero-field cooled to the lowest T: (a) Sam- 
ples 1 and 6; (b) samples 2, 4, 4A and 4B; (c) samples 3, 5 
and 7. 

X(T)H Y 2 - y 3 (y w 0.03), on the other hand, sus- 
tains its high-T metallic state down to low tempera- 
tures, and at its Neel temperature Tn ~ 10 K it under- 
goes a transition to an antiferromagnetic phase with a 
cusp in x(T).l3 V3O5 also orders antiferromagnetically 
at Tn = 75.5 K, but_^(T) shows a broad maximum at 
a higher T = 125 Ko Though not detected in our x- 
ray diffraction measurements, V4O7, which has the same 



V oxidation state as in LiV204, also displays a cusp in 
X ohs (T) at T N ps 33 K and x° bs (T) follows the Curie- 
Weiss law for T > 50K.E-3 The susceptibilities of these 
V-0 phasps.|ajfe all on the order of 10~ 4 to 10 -3 cm 3 /mol 
at low TEM Moreover, the T variations of x° hs (T) in 
these vanadium oxides for T < 10 K are, upon decreasing 
T, decreasing (V2-VO3) or nearly T independent (V3O5 
and V4O7), in contrast to the increasing behavior of our 
Curie-like impurity susceptibilities. From the above dis- 
cussion and the very small amounts of V-0 impurity 
phases found from the Rietveld refinements of our x-ray 
diffraction measurements, we conclude that the V-0 im- 
purity phases cannot give rise to the observed Curie-like 
upturns in our x° hs {T) data at low T. These Curie-like 
terms therefore most likely arise from paramagnetic de- 
fects in the spinel phase and/or from a very small con- 
centration of an unobserved impurity phase. 

Figure ^(b) shows how the additional heat treatments 
of the as-prepared sample 2 yield different behaviors of 
X obs (T) at low T in samples 4, 4A and 4B. Only liquid- 
nitrogen quenching (sample 4) caused a decrease in the 
Curie- like upturn of sample 2. On the contrary, ice wa- 
ter quenching (sample 4A) and oven-slow cooling (sam- 
ple 4B) caused x° bs (T) to have an even larger upturn. 
However, the size of the Curie-like upturn in x° hs (T) of 
sample 4 was found to be irreproducible when the same 
liquid-nitrogen quenching procedure was applied to an- 
other piece from sample 2; in this case the Curie-like 
upturn was larger, not smaller, than in sample 2. The 
observed susceptibility (not shown) of this latter liquid 
nitrogen-quenched sample is very similar to those of sam- 
ples 4A and 4B. The x° bs (T) of pamnles 4 A and 4B re- 
semble those reported previouslyJi3c3 



2. Isothermal Magnetization versus Magnetic Field 

Larger Curie-like upturns were found in samples with 
larger curvatures in the isothermal M ohs (H) data at low 
T. A few representative M ohs (H, 2K) data for samples 
showing various extents of curvatures in M° hs (H) are 
shown in Fig. El which may be compared with the corre- 
sponding x° hs \T) data at low T in Figs. |[ This correla- 
tion suggests that the Curie- like upturns in x° bs (T) arise 
from paramagnetic (field-saturable) impurities/defects in 
the samples. On the other hand, there is no obvious 
correlation between the magnetic impurity concentration 
and the V2O3 or V3O5 phase impurity concentration, as 
noted above. 

The isothermal M ohs {H) data for H < 5.5 T displayed 
negative curvature for T < 10-20 K and linear behavior 
for higher T, as illustrated for sample 1 in Fig. [|. The 
concentrations and other parameters of the magnetic im- 
purities in the various samples were obtained from anal- 
yses of M° bs (H) isotherms as follows. From high-field 
measurements, the intrinsic magnetization M(H, 0.5 Kl 
of L1V2O4 is proportional to H up to H ~ 16T.CJ 
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Therefore, the observed molar mag netization M ohs (H,T) 
isotherm data for each sample were fitted by the equation 

M° hs (H, T) = M imp (H, T) + M(H, T) 

= fim P N A gi m ptJ.BSi m pBs ilnp (x) +x(T)H , 

(1) 
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FIG. 4. Comparison of the negative curvatures of observed 
magnetization isotherms M° hB at T = 2 K versus applied mag- 
netic field H for samples 1, 3, 5 and 7. 




FIG. 5. Observed magnetization M° bs versus applied mag- 
netic field H isotherms at temperatures T = 2, 10 and 300 K 
for LiV204 sample 1. Negative curvature in M ohB (H) is not 
present for T > 10 K for this sample. 

where /imp is the magnetic impurity concentration, Na 
Avogadro's number, <?i mp the impurity g-factor, /ib 
the Bohr magneton, S lmp the impurity spin, i?s imp 
the Brillouin function, x the intrinsic susceptibility of 
the L1V2O4 spinel phase and H the applied magnetic 
field. The argument of the Brillouin function is x = 
gim V mSi mp H /[k B (T - 6> imp )]. 6> imp represents the Weiss 
temperature of the Curie- Weiss law when the suscepti- 
bility is obtained by expanding the Brillouin function in 
the limit of small H/(T — #i mp ). Incorporating the pa- 
rameter 9 lmp ^ takes account of possible interactions 



between magnetic impurities in a mean-field manner. To 
improve the precision of the obtained fitting parameters, 
we fitted M° hs (H) isotherm data measured at more than 
one low temperature simultaneously by Eq. (^). Since the 
negative curvature of the isothermal M ohs (H, T) data di- 
minishes rapidly with increasing T, only low T (1.8-6 K) 
data were used. Furthermore, a linear T dependence of 
x(T) in this T range was assumed [see Fig. |3|(a)] in or- 
der to reduce the number of free parameters. However, 
x(T = 2K) and the linear slope dx/dT still have to be 
determined. Hence up to six free parameters were to be 
determined by fitting Eq. (|l|) to the data: /imp, ffimp, 



Simp 7 ^imp j 



X (T = 2K) and dx/dT. 



With all six parameters varied as free parameters, fits 
of M ohs (H, T) by Eq. (|l|) produced unsatisfactory results, 
yielding parameters with very large estimated standard 
deviations. Therefore, we fixed Sj mp to various half- 
integer values starting from 1/2, thereby reducing the 
number of free parameters of each fit to five. With re- 
gard to the gimp values, g-factors of slightly less than 2 
are observed in V +4 compounds: VO2 (1.964) (Ref. [l6L 
(NH 4 ) :E V 2 05 (1.962) (Ref. ^) and U x Y 2 O b (1.96)13 
Using gimp ~ 2 as a guide, we selected a few values 
of >Simp which resulted in g ~ 2 in the five-parameter 
fit. Then using the obtained parameter values we cal- 
culated and plotted the impurity magnetization Mi mp 
(= M obs - X H) versus H/(T - 9 imp ) for all the low T 
data utilized in the fit by Eq. (Q). If a fit is valid, then 
all the Mi mp [H/(T — 6*imp)] data points obtained at the 
various isothermal temperatures for each sample should 
collapse onto a universal curve described by Mi mp = 
fimpN A gim P tJ>BSi mp Bs iltlp (x). The fixed value of SUp 
which gave the best universal behavior for a given sample 
was chosen. Then, using this Simp, we fixed the value of 
.gimp to 2 to see if the resultant M- lmp [H/{T — 6* im p)] data 
yielded a similar universal behavior. For the purpose 
of reducing the number of free parameters as much as 
possible, if this fixed-g fit did yield a comparable result, 
the parameters obtained were taken as the final fitting 
parameters and are reported in this paper. For sam- 
ple 1 only, the fit parameters obtained by further fixing 
#im P = are reported here. To estimate the goodness of 
a fit, the x 2 P er degree of freedom (DOF) was obtained, 



which is defined as (JV P - P)" 1 " ^ alc ) 2 /^ 2 , 

where N p is the number of data points, P is the number 
of free parameters, and Oi is the standard deviation of 
the observed value Mj. A fit is regarded as satisfactory 
if % 2 /DOF < 1, and this criterion was achieved for each 
of the nine samples. 

The magnetic parameters for each sam ple, obtained as 
described above, are listed in Table III. Plots of Mi mp 
versus H/(T — #i mp ) for the nine samples are given in 
Figs. ||(a), (b) and (c), where an excellent universal be- 
havior for each sample at different temperatures is seen. 
The two magnetically purest samples 1 and 6 have the 
largest relative deviations of the data from the respective 
fit curves, especially at the larger values of H / (T — 6>i mp ). 
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Since these two samples contain extremely small amounts 
of paramagnetic saturable impurities, the magnetic pa- 
rameters of the impurities could not be determined to 
high precision. The impurity spins i5i mp obtained for 
the nine samples vary from 3/2 to 4. In general, the 
magnetic impurity Weiss temperature |6*i mp | increased 
with magnetic impurity concentration /imp- From the 
chemi- cal analyses of the starting materials (V2O5, 
NH4VO3 and Li 2 C03) supplied by the manufacturer, 
magnetic impurity concentrations of 0.0024 mol % Cr and 
0.0033 mol%Fe are inferred with respect to a mole of 
L1V2O4, which are too small to account for the param- 
agnetic impurity concentrations we derived for our sam- 
ples. 



3. Magnetization versus Temperature Measurements 

a. Low Magnetic Field ZFC and FC Measurements 
The zero-field-cooled (ZFC) x° hs (T) data in Fig. |(a) for 
our highest magnetic purity samples 1 and 6 show a broad 
maximum at T peak « 16 K. One interpretation might be 
that static short-range (spin-glass) ordering sets in be- 
low this temperature. To check for spin-glass ordering, 
we carried out low-field (10-100 G) ZFC and field-cooled 
(FC) magnetization measurements from 1.8-2 K to 50 K 
on all samples except samples 2 and 4B. For each sample, 
there was no hysteresis between the ZFC and FC mea- 
surements, as illustrated for sample 4 in Fig. |7[ ancLthus 
no evidence fopspin-glass ordering above 1.8-2 K.C3 

Ueda et a/.E3 reported that spin-glass ordering oc- 
curs in the zinc-doped lithium vanadium oxide spinel 
Lii_ 2 .Zn :c V204 for 0.1 < x < 0.9. However, spin-glass 
ordering was not seen in the pure compound L1V2O4, 
consistent with our results. Further, positive-muon spin 
relaxation /iSR measurements for sample 1 did not de- 
tect static magnetic ordering down to 20 mK.El However, 
the /xSR measurements did indicate the presence of static 
spin-glass^prdering in the off-stoichiometric sample 3 be- 
low 0.8 K.Q As mentioned in Sec. Ill A, the stoichiometry 



of sample 3 was intentionally made slightly cation- 



TABLE III. Results of magnetization M° hs (H,T) isotherm analyses, where the T values used are listed in the second 
column, /mag imp is the molar magnetic impurity concentration. The error in the last digit of a quantity is given in paren- 
theses. All numbers without an error listed were fixed in the fit. The Curie constant of the impurities was calculated from 

Cimp = /"mag imp^VAf/imp/^B 5imp (5*imp + l)/(3fc B ). 
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FIG. 7. Observed magnetic susceptibility x° bs ( T ) = 
M° hs {T)/H versus temperature T in a low magnetic field 
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and in the low field (FC). 
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6 and (b) samples 2, 3, 4, 5, 7, 4A and 4B. The solid lines are 
guides to the eye. 

deficient, and may contain cation vacancies. Such a de- 



fective structure could facilitate the occurrence of the 
spin-glass behavior by relieving the geometric frustration 
among the V spins. Whether the nature of the spin-glass 
ordering in sample 3 is similar to or different from that 
in Lii_ a; Zn 2 ;V204 noted above is at present unclear. 

b. Intrinsic Susceptibility The intrinsic susceptibil- 
ity x{T) was derived from the observed M ohs (T) data at 
fixed H — IT using X (T) = [M ohs (T) - M imp (H 7 T)]/H, 



where M imp (H, T) is given by Eq. (m) with H — IT and 



by the parameters for each sample given in Table III , and 
T is the only variable. The x(T) for each of the nine sam- 
ples is shown in Figs. ||(a) and (b), along with x° bs (D 
for samples 1 and 6. A shallow broad peak in x(T) is 
seen at a temperature T p0 ak = 18, 16, 18, 18, 15, 17, 17, 5 
and 14 K for samples 1-7, 4A and 4B, respectively. The 
peak profiles seen in x(T) for the two magnetically purest 
samples 1 and 6 are regarded as most closely reflecting 
the intrinsic susceptibility of L1V2O4. This peak shape 
is obtained in the derived X (T) of all the samples ex- 
cept for sample 4 A, as seen in Fig. ^(b). The physical 
nature of the magnetic impurities in sample 4A is evi- 
dently different from that in the other samples. Except 
for the anomalous sample 4A, the X (T = 0) values were 
estimated from Figs. |8|(a) and (b), neglecting the small 
residual increases at the lowest T for samples 2, 6, 7 and 
4B, to be 

X(0)=9.8, 10.8, 9.6, 9.7, 10.0, 10.2, 10.2, 

9.8xl0~ 3 cm 3 /mol (samples 1-7, 4B) . (2) 



IV. MODELING OF THE INTRINSIC 
MAGNETIC SUSCEPTIBILITY 

A. The Van Vleck Susceptibility 

The Van Vleck paramagnetic orbital susceptibility x vv 
may be obtained in favorable cases from the so-called K- 
X analysis, i.e., if the transition metal NMR frequency 
shift K depends linearly on Xi with T an implicit pa- 
rameter. One decomposes x(^) P er mole of transition 
metal atoms according to X (T) = x COIC + X VV + X spin (T). 
We neglect the diamagnetic orbital Landau susceptibil- 
ity, which should be small for (i-electron bands Ei The 
NMR shift is written in an analogous fashion as 

K{T) = K VY + K spin (T) ; (3) 

a term K cmc does not appear on the right-hand side of 
Eq. (|^) because the absolute shift due to ^ corc is expected 
to be about the same as in the Knight shift reference com- 
pound and hence does not appear in the shift measured 
with respect to the reference compound. Each compo- 
nent of K is written as a product of the corresponding 
component of x and of the hyperfinc coupling constant 
A as 
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K vv = 

^spin _ 



NaMb' 

y/jspin 



,W 



-X 



Combining Eqs. (||) and (||) yields 

/(VV A spin 

X = _ x vv + x s p- 



(4a) 
(4b) 

(5) 



If K(T) varies linearly with x(T), then the slope is 
A spm / Na/J-b since x VV (and x COIC ) is normally indepen- 
dent of T. We write the observed linear relation as 



/(spin 

K = K + X 

N a /J>b 



Setting the right-hand-sides of Eqs. (||) and (||) equal to 
each other gives 



,vv 



N A y B K + A s P in x c 
A vv - A s P in 
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FIG. 9. 51 V NMR Knight shift K versus obseryed-Siagnetic 
susceptibility )fn s for L1V2O4 by Amako et ai£3'E3 and by 
Mahajan et alxB for LiV2C>4 sample 2. The lines are linear 
fits to the data according to Eq. (H). 



From 51 V NMR and x(T) measurements, the K vs. 
X relationship for LiV204 was determined by Amako et 
al.Ln and was found to be linear from 100-300 K, as shown 
in Fig. |9|. Our fit to their data gave 



K = 0.0117(4) 



17.08(21) 



mol V 



/ cm 3 



mol V 



(8) 



shown as the straight line in Fig. |^. Comparison of 
Eqs. (§) and (§) yields 



#0 = 0.0117(4) , 
A spin = -95.4(12) kG 



(9a) 
(9b) 



The orbital Van Vleck hyperfine coupling constants 
for V +3 and V +4 are similar. For atomic V +3 , one has 
A vv = 403 kG (Ref. |lL We will assume that A vv in 
LiV204 is given by that£3 for atomic V +4 , 



A vv = 455kG 



(10) 



The core susceptibility is estimated here from Sel- 
wood's table,E3 using the contributions [in units of 
-10- 6 cm 3 /(mol ion)] 1 for Li+\ 7 for V+ 4 and 12 for 
0~ 2 , to be 



= -63 x 10" 



cm 
mol 



(6) Inserting Eqs. ©-© into (g) yields 



X 



vv 



2.48(9) x 10~ 4 ^- 
mol 



(11) 



(12) 



Mahajan et aim have measured the 51 V K(T) for our 
LiV204 sample 2 from 78 to 575 K. Their data are plot- 
ted versus our measurement of x (T) for sample 2 from 
74 to 400 K in Fig. l| Applying the same K-\ analysis 
as above, we obtain 



K = 0.0101(3) , 
A spm = _ 76 .9( 8 ) kG 



X vv = 2.22(6) x 10- 4 ^- 
mol 



(13) 
(14) 

(15) 



where the linear fit of K vs. x° hs is shown by the dashed 
line in Fig. ^. 

We may compare our similar values of x VV for LiV204 
in Eqs. ( |l2|) and ([l5|) with those obtained from K- 
X analyses of other oxides containing V +3 and V +4 . 
For stoichiometric V2O3 above its metal-insulator tran- 
sition—temperature of ~ 160 K, JonesEfJ and Takigawa 
et alxA respectively obtained x VV i-t= 2.10 and 2.01 x 
10- 4 cm 3 /(mol V). Kikuchi et alS3 obtained x VV = 
0.92 *-40~ 4 cm 3 /(mol V) for LaV0 3 , and for V0 2 , Pouget 
et alE3 obtained x VV = 0.65 x 10- 4 cm 3 /(mol V). 



B. High- Temperature Series Expansion Analysis of 
the Susceptibility 

Above ~ 50 K the monotonically decreasing suscepti- 
bility of LiV204 with increasing T has been interpreted 
by previous workers in terms of the Cucie-|Weiss law for 
a system of spins S = 1/2 and g « 2.t2rii3 To extend 
this line of analysis, we have fitted x(T) by th^rbigh- 
temperature series expansion (HTSE) prediction! 56 !! 57 ! up 
to sixth order in l/T. The assumed nearest-neighbor 
(NN) Heisenberg Hamiltonian between localized mo- 
ments reads TL = JJ2(ij) Si'Sji where the sum is over all 



10 



NN pairs, J is the NN exchange coupling constant and 
J > denotes AF interactions. A HTSE of x ^ n SE (T) 
arising from this Hamiltonian up to the n max -th order of 
J/fc B T for general lattices and spin S was determined by 
Rushbrooke and Wood,K3 given per mole of spins by 



3fc B T 



( J 



(16) 



where 60 = 1- The b n coefhcients for S = 1/2 up to sixth 
order (n max = 6) are 

Z Z Z ^Pl_ 

h = b 2 = -, 6 3 = -(l-— ) , 



15360 



184320 



15pi + 5p 2 ) , 
(90z - 122 + 245pi - 60z Pl - 45pj 

-90p 2 + 25p 3 ) , 
(134z 2 - 783z + 713 + 908z Pl - 2697pi 

- I06zpl + I284pj - 234zp 2 

+ 849p 2 - 291p 3 + 75p 4 - 288pip 2 

-51g-8r) . (17) 



Here z is the nearest-neighbor coordination number, and 
p n , q and r are so-called lattice parameters which de- 
pend upon the geometry of the magnetic lattice. The 
Curie Law corresponds to maximum order n max = 0, and 
the Curie- Weiss Law to maximum order n max = 1. For 
the B sublattice of a normal-spinel structure compound 
yl[£> 2 ]04, which is geometrically frustrated for AF inter- 
actions, the parameters are z = 6, p\ = 2, p 2 = 2, p% — 0, 
Pi = 2, q = 0, r = 2. For S = 1/2, Eq. @ then yields 



h = 



3 

2 ' 

_ 37 
~128 



b 2 = 



3 

4 ' 
h = 



bn = - 



1 

16 ' 



43 
640 



6 fi = 



1361 
6144 



(18) 



Figure |To| illustrates the HTSE predictions of Eq. ( |l6|) for 
S = 1/2 using these 6 n coefficients for n max = 1 to 6. The 
theoretical Xhtse(^) predictions with n max — 2, 3 and 6 
exhibit broad maxima as seen in our experimental x(T) 
data. The prediction with n max — 6 is evidently accu- 
rate at least for fc B X '/ 'J > 1.6; at lower T, the theoretical 
curves with 7-i max = 5 and 6 diverge noticably from each 
other on the scale of Fig. |lC]. Our fits given below of the 
experimental data by the theoretical Xhtse(-^) predic- 
tion were therefore carried out over temperature ranges 
for which fc B T/ J > 1.6. The Weiss temperature 9 in the 
Curie- Weiss law is given for coordination number z = 6 
and S = 1/2 by 9 = -zJS{S + l)/(3fc B ) = -3J/(2fc B ). 

To fit the HTSE calculations of Xhtse(^) to experi- 
mental data, we assume that the experimentally deter- 
mined intrinsic susceptibility x(T) is the sum of a T- 
independent term \o and Xhtse(-^)' 
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FIG. 10. High-temperature series expansion predictions 
of the normalized spin susceptibility Xhtse'-V-^A-S' 2 A. with 
11™ = 1-6 versus reduced temperature k B T/J [Eq. @] 
for the antiferromagnetically coupled spins S = 1/2 in the B 
sublattice of a normal-spinel compound A[B2]C>4. 
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FIG. 11. Intrinsic susceptibility \ versus temperature T for 
LiV2 04 sample 1 (filled circles) and fits (curves) by the high-T 
series expansion (HTSE) prediction to 6th order in 1/T for 
the 50-400 and 100-400 K temperature ranges. 



X(T) = Xo + Xhtse( T ) 



(19) 



with Xhtse(-^') given by Eq. ( |l6| ) and the b n coefficients 
for S = 1/2 in Eq. (|l8|). The three parameters to be 
determined are xoj 9 and J/fc B - The fitting parameters 
for samples 1-7, 4A and 4B using n max = 6, and for 
sample 1 also using ^ max — 2 and 3, are given in Table IV 
for the 50-400 and 100-400 K fitting ranges. The fits for 
these two fitting ranges for sample 1 and n max = 6 are 
shown in Fig. [ll|. Both g and J/fc B tend to decrease 
as the lower limit of the fitting range increases. The 
HTSE fits for all the samples yielded the ranges C — 
N A g 2 ^y (4k B ) = 0.36-0.48 cm 3 K/(molV) and 9 = -20 
to — 42 K, in agreement with those reported previously 
(see Table |) . xo was found to be sensitive to the choice 
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of fitting temperature range. For the 50-400 K range, 
Xo was negative for some samples. Recalling the small 
negative value of the core diamagnetic contribution in 
Eq. ([ll]) and the larger positive value of the Van Vleck 
susceptibility in Eqs. (O) and (llq), it is unlikely that xo 
[defined below in Eq. fl2CQl would be negative. Negative 
values of xo occur when the low-T limit of the fitting 
range is below 100 K, and may therefore be an artifact of 
the crossover between the local moment behavior at high 
T and the HF behavior at low T. 

To eliminate xo as a fitting parameter, we also fitted 
dx/dT by the HTSE prediction for that quantity. The 
experimental dx/dT was determined from a Pade ap- 
proximant fit to x(T) an d is plotted in Fig. [l2] for sam- 
ple 1. These data were fitted by dx^fxsE/^ obtained 
from the HTSE prediction Eq. (|l^) with n max = 6, where 
the fitting parameters are now g and J/fce- The fits were 
carried out over the same two T ranges as in Fig. [Tl| ; 
Table ^ displays the fitting parameters and the fits are 
plotted in Fig. [l2[ Both g and J/fce were found to be 
larger than the corresponding values in Table IV . Of the 



two fitting ranges, the 100-400 K fit is the best fit inside 
the respective range, though it shows a large deviation 
from the data below this range. Using the fitting pa- 
rameters, the HTSE x spin (T) is obtained from Eq. @. 
According to Eq. (p9|), the difference between the exper- 
imental x(T) and Xhtse( t )> S x(T) = x{T) - Xhtse( t X 
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FIG. 12. Temperature derivative of the experimental in- 
trinsic susceptibility, dx/dT, for LiV204 sample 1 (heavy 
solid curve). Fits by the T derivative of the HTSE predic- 
tion dx Bpln /dT in Eq. ( |l6| ) are also shown for T ranges of 
50-400 K (dashed curve) and 100-400 K (light solid curve). 

should represent the T-independent contribution xo- 
Sx(T) is plotted for sample 1 versus T in Fig. [p| for 
the 50-400 K and 100-400 K fit ranges. Again, the su- 
periority of the 100-400 K fitting range to the other is 
evident, i.e., xo is more nearly constant for this fitting 
range, xo f° r the 50-400 K fit range is negative within 
the range. This sign is opposite to that obtained in the 



TABLE IV. Results of high-temperature series expansion calculation fits to the intrinsic magnetic susceptibility data for 
LiV204 over the temperature ranges 50-400 K and 100-400 K. The error in the last digit of a quantity is given in parentheses. 



Sample 


n max 




50-400 K 






100-400 K 




No. 




Xo 


g 


J/k B 


Xo 


9 


J/k B 






(10~ 4 cm 3 /mol) 




(K) 


(10~ 4 cm 3 /mol) 




(K) 


1 


2 


0.8(4) 


2.17(1) 


25.8(5) 


2.7(3) 


2.07(2) 


20(1) 


1 


3 


0.7(4) 


2.18(2) 


26.2(6) 


2.6(3) 


2.07(2) 


20(1) 


1 


6 


0.5(4) 


2.19(2) 


26.9(7) 


2.6(3) 


2.07(2) 


20(1) 


2 


6 


-0.2(5) 


2.26(2) 


26.7(8) 


2.6(3) 


2.11(2) 


19(1) 


3 


6 


-1.3(5) 


2.23(2) 


27.8(7) 


1.4(3) 


2.08(2) 


20.5(8) 


4 


6 


1.1(6) 


2.16(3) 


26.4(9) 


4.1(5) 


1.99(3) 


17(2) 


4A 


6 


-0.6(8) 


2.20(3) 


26(1) 


2.3(2) 


2.05(1) 


18.1(6) 


4B 


6 


-0.7(5) 


2.12(2) 


26.2(8) 


1.8(5) 


1.97(3) 


18(2) 


5 


6 


1.2(7) 


2.17(3) 


25(1) 


4.9(7) 


1.95(4) 


13(2) 


6 


6 


0.8(1) 


2.251(6) 


26.5(2) 


3.3(7) 


2.108(4) 


18.4(2) 


7 


6 


0.5(3) 


2.20(1) 


25.8(5) 


3.0(1) 


2.051(8) 


17.5(4) 



TABLE V. Parameters g and J/k B obtained by fitting the temperature T derivative of the experimental intrinsic suscepti- 
bility data for LiV204 samples 1 and 6 by the T derivative of the HTSE spin suscepti bility [Eq. (0)] with n max = 6 for two 
different temperature ranges of the fit. The T-independent susceptibility xo was determined by averaging Sx(T), see Fig. [p| 
The error in the last digit of a quantity is given in parentheses. 







50-400 K 






100-400 K 




Sample No. 


Xo 


9 


(J/k B ) 


Xo 


9 


{J/ks) 




(io~ 4 f£) 




(K) 


(HT 4 £si) 

mol / 




(K) 


1 


-1.5(1) 


2.275(3) 


29.61(7) 


2.00(4) 


2.103(2) 


22.27(8) 


6 


-2.73(5) 


2.402(4) 


31.61(9) 


2.11(1) 


2.174(3) 


22.1(1) 
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first HTSE fitting results in Tabic [V . This inconsistency 
found in the fit using a low T limit below 100 K may 
again be due to changing physics in the crossover regime, 
which would invalidate the parameters. By averaging the 
Xo values for samples 1 and 6 in the given ranges, we 
obtained the T-independent contribution xoj as listed in 
Table [v|. 
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FIG. 13. The differences between the experimental intrin- 
sic susceptibility x(T) of L1V2O4 sample 1 and the HTSE 
prediction x spm obtained from the T derivative analysis, 
5x{T) = x(T) - X sph \ versus temperature T for the fit- 
ting T ranges of 50-400 K (open squares) and 100-400 K 
(filled circles). For a valid fit, these differences should be 
the T-independent susceptibility \o- 

In the itinerant plus localized moment model implicitly 
assumed in this section, xo can be decomposed as 



vv 



Xo = X^ + X vv +X 



Pauli 



(20) 



where ^ Paul1 is the temperature-independent Pauli spin 
susceptibility of the conduction electrons. Using the re- 
sults of x corc [Eq. ©], x VV [Eq. ©] and Xo (100-400 K 
range, Table M), we find 



Pauli 



Pauli 



0.41(10) x 10~ 4 cm 3 /mol (sample 1) , (21a) 



0.52(7) x 10~ 4 cm 3 /mol (sample 6) . (21b) 



These ^- Pauh values are approximately four timeSj-smallcr 
than that obtained for LiV204 by Mahajan et altB They 
used x° hs (T) in the T range 100-800 K, combining o 
X obs (T) data to 400 K with those of Hayakawa et al. 
to 800 K. By fitting these combined data by the expres- 
3hs (T) = X o + 2C/(T-6), they obtained 



Xo = 



sion x" 

5.45 x 10 -4 cm 3 /mol. As shown above and also discussed 
in Ref. (3f| the value of xo is sensitive to the fitting tem- 
perature range. For LiTi 2 04, ^ Pauh ~2x 10~ 4 cm 3 /mol 
(Refs. |23j , |58| ) between 20 and 300 K, which is a few times 
larger than we find for L1V2O4 from the 100-400 K range 
fits (Table 0). 



C. Crystal Field Model 

The ground state of a free ion with one 3d electron 
is 2 D 3 / 2 and has five-fold orbital degeneracy. The point 
symmetry of a V atom in L1V2O4 is trigonal. If we con- 
sider the crystalline electric field (CEF) seen by a V atom 
arising from only the six nearest-neighbor oxygen ions, 
the CEF due to a perfect oxygen octahedron is cubic (Oh 
symmetry), assuming here point charges for the oxygen 
ions. In this CEF the degeneracy of the five d orbitals 
of the vanadium atom is lifted and the orbitals are split 
by an energy "lODg" into a lower orbital t 2g triplet and 
a higher orbital e g doublet. However, in LiV204 each 
V-centered oxygen octahedron is slightly distorted along 
one of the < 1 1 1 > directions [see Fig. 0(b)], as discussed 
in Sec. Ill A. This distortion lowers the local symmetry 
of the V atom to (trigonal) and causes a splitting of 
the tig triplet into an A\ g singlet and an E g doublet. It 
is not clear to us which of the E g or A\ g levels become 
the ground state, and how large the splitting between the 
two levels is. These questions cannot be answered readily 
without ^-knowledge of the magnitudes of certaia, radial 
integrals,t3 and are not further discussed here.E2l How- 
ever, this trigonal splitting is topically about an order of 
magnitude smaller than lQDqSA In the following, we will 
examine the predictions for x{T) of a d 1 or d 2 ion in a 
cubic CEF and compare with our experimental data for 

LiV2 ° 4 - B 

KotaniErl calculated the effective magnetic moment 
fi e g = PcffMB per G?-atom for a cubic CEF using the Van 
Vleck formula£3 The spin-orbit interaction is included, 
where the coupling constant is A. For an isolated atom 
/ioff(T) is defined by x(T) = N^ e {T)/(3k B T), where 
/U e ff is in general temperature-dependent and N is the 
number of magnetic atoms. With spin included, one uses 
the double group for proper representations of the atomic 
wavefunctions. Then in this cubic double group with one 
d-electron the six-fold (with spin) degenerate t^-.lejfel 
splits into a quartet T s (t 2g ) and a doublet r 7 (i 2ff )Be5El 
The four-fold degenerate e g level does not split and its 
representation is Ts{e g ). For a positive A, as is appro- 
priate for a 3d atom with a less than half-filled d-shell, 
^8{t2 g ) is the ground state, and the first-order Zeeman ef- 
fect does not split it; this ground state is non-magnetic. 
Kotani does not include in his calculations of /x e ff the 
possible coupling of Fg(t2 g ) and rg(e 9 ), which have the 
same symmetry, and assumes that the cubic CEF split- 
ting lODq is large enough to prevent significant mix- 
ing. On the other hand, the cubic double group with 
two d-electrons gives an orbitally nondegenerate, five- 
fold spin-degenerate, ground state with angular momen- 
tum quantum number J = 2 which splits into five non- 
degenerate levels under a magnetic field. The spin-orbit 
coupling constant is A =^-J- 250 cm -1 for d 1 (V +4 ) and 
+105 cm -1 for d 2 (V +3 ).El The effective moment is de- 
fined from the observed molar susceptibility of LiV2C>4 as 
X° hs (T) - Xo + 2N A [p°* s (T)} 2 f i 2 B /{3k B T), where we take 
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Xo = 2.00 x 10 4 cm 3 /mol given in Table [v|. Kotani's 



results from the Van Vleck equations are 



Pctt 



8+ (3a: - 8)e" 



x{2- 



1/2 



(22) 



for the d 1 ion, and 



PcS 



3(§x + 15 + (f + 9)e- x - 24e-% x ) 



"I 1/2 



x(5 + 3e 



(23) 



for the d 2 ion, where x = A/fesT. Figure |lj shows p°g s , 
and as a function of T. For comparison is also 

shown _Pgg +2 ' ) obtained by assuming that p°g s (T) arises 
from an equal mixture of V +3 and V +4 localized mo- 
ments. None of the three calculated curves agree with 
the experimental data over the full temperature range. 
However, in all three calculations p c ff increases with T, 
in qualitative agreement with the data, perhaps implying 
the importance of orbital degeneracy in L1V2O4 and/or 
antifcrromagnetic coupling between vanadium spins. The 
nearly T-independent p°% s w 1.8 for T > 100 K is close to 
the spin-only value p e s = g^J S(S + 1) with 5 = 1/2 and 
g m 2, as expected in the absence of orbital degeneracy; 
however, this result also arises in the theory for the d 1 
ion when k^T ~ A, as seen by comparison of the solid 
curve with the data in Fig. Il4| at ~ 300 K. 




150 200 250 300 350 400 
T(K) 

FIG. 14. Observed effective magnetic moment in fts, P°ff S , 
versus temperature T of LiV2 04 sample 1 (filled diamonds). 
Also shown as the curves are the predictions for d 1 ions 
and £Tg for d 2 ions by Kotani,E3 and pi for an equal mix- 
ture of d 1 and d 2 ions, in a cubic crystalline electric field, 
including spin-orbit coupling. 



D. Spin-1/2 Kondo Model and Coqblin-Schrieffer 
Model 

X (T) data for /-electron HF compounds are often 
found to be similar to the predictions of the single-ion 



jantion 
3 The 



Kondo modelESB for spin S =1/2 or its 
to S > 1/2 in the Coqblin-Shrieffer model 
zero-field impurity susceptibility xcs (T) of the Coqblin- 
Shrieffer model was calculated exactly as a function of 
temperature by RajanO His numerical results xcs(T) 
for impurity angular momentum quantum number J = 
1/2, 7/2 show a Curie- Weiss-like 1/T dependence 
(with logarithmic corrections) for T 3> Tk, where Tk 
is the Kondo temperature. As T decreases, Xcs(T) 
starts to deviate from the 1 /T dependence, shows a peak 
(at T w 0.2T K ) only for J > 3/2, and levels off for 
T < 0.2T K for all J. 

In the zero temperature limit the molar susceptibility 
for J = S = 1/2 (which corresponds to the S = 1/2 
Kondo model) isE3 



Xcs(T = 0) 



0.102678A r A5>B 
kv,T K 



(24) 



Setting g = 2, and using the intrinsic x(T ~ > 0) 
0.0049 cm 3 /(molV) for LiV 2 4 sample 1 from Eq. ([ 
Eq. (^J) yields the Kondo temperature 



T K = 32.1 K 



(25) 



On the other hand, if the g-value of 2.10 from Table |v| 
(100-400 K range) is employed instead, the Kondo tem- 
perature is 



T K = 35.5 K 



(26) 



The temperature dependence of the impurity suscepti- 
bility of the S = 1/2 Kondo model was obtained using ac-. 
curate Bethe ansatz calculations by Jerez and Andrei.O 
Their T — > value for the coefficient on the right-hand- 
side of Eq. (||) is 0.1028164, about 0.1% too high com- 
pared with the correct prefactor in Eq. (|2~i|). We fitted 
their calculated values for t = 0.00104 to 102.53 by 



4xcsfcBr 
Ng^l 



1 + T + w + # 



4(0.1028164)ti 5 
t 5 



1 + f 



d 2 

7^ 



— 



(27a) 



ni = 530.417 , 7i 2 = 4697.91 , n 3 = 1404.18 , 
n 5 = -418.781 , di = 695.557 , d 2 = 8605.97 , 
d 3 = 11373.7 , di = 2937.88 , (27b) 



where t = T/Tk- Equation (27a) has the correct form 
Xcs(0) + bt 2 at low T and approaches a Curie law in 
the high-T limit, as required by the Kondo model. The 
large rij and di coefficients arise because Xcs(T) con- 
verges very slowly to the Curie law at high tempera- 
tures. The rms deviation of the fit values from the Bethe 
ansatz calculation values is 0.038%, and the maximum 
deviation is 0.19% at t = 66.9. Using the above-stated 
g- values and T K from Eqs. @ and @, the S = 1/2 
Yes (T) calculations are c omp ared with our x(T) data in 
Fig. [15]. Note that in Fig. 15, both the T-independent xo 
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(Table |yj) and impurity susceptibilities are already sub- 
tracted from x obs - Although the T K values in Eqs. 
and (|2^) are comparable to those obtained from specific 
heat analyses, Bo the S = 1/2 Kondo model predictions 
for x(T) with these 7k values do not agree with our ob- 
served temperature dependence. This failure is partly 
due to the fact that our x(T) data exhibit a weak maxi- 
mum whereas the S = 1/2 Kondo model calculation does 
not. 
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LiV 2 4 Sample 1 



y [g=2, T =32.1 K] 
X [q=2.103, T =35. 5K] 




50 100 150 200 250 300 350 400 
T(K) 

FIG. 15. Temperature T-dependent part of the mag- 
netic susceptibility, x ~ Xo, versus T for LiV2 04 sample 
1 (filled circles). Also shown as solid and dashed curves 
are the predictions of the spin 5=1/2 Kondo model for 
(g, Ik) = (2, 32. IK) and (2.103, 35. 5K), respectively, where 
g is the p-factor and 7k is the Kondo temperature. 



As noted above, the Coqblin-J 
3/2 does give a peak in xcs(7').l 



shrieffer model for J > 
i Defining the ratio 



r(%) = 100 



peak 

Xcs 



Xcs(0) 



Xcs(0) 



(28) 



where 



peak 



calculationsO give r = 2, 7, 



is the value of Xcs(T) at the peak, the 
11, 17 and 22% for J = 3/2, 
2, 5/2, 3 and 7/2, respectively. The observed value is 
r = 8.2% in sample 1, which is between the theoretical 
values for J = 2 and 5/2. Fits of xcs{T) to our x(T) 
data of sample 1 for T = 2-400 K are shown in Fig. [L6| 
and the parameters are 

Xo = 2.3(3) x 10" 4 cm 3 /mol , 
.9 = 0.790(3), 7k = 97.8(6) K (J = 2) ; (29a) 

Xo = 6.9(9) x 10" 4 cm 3 /mol , 
g = 0.591(7) , 7k = 103(2) K (J = 5/2) . (29b) 

The J — 2 curve fits our x(T) data fairly well. How- 
ever, the 1.5 d-electrons per V ion could not give rise to 
a J value this large; the very small value of g is also 
considered highly unlikely. 

On the basis of the above analysis we conclude that the 
Coqblin-Schrieffer model for S > 1/2 and the 5=1/2 



Kondo model cannot explain the intrinsic susceptibility 
of L1V2O4 over any appreciable temperature range. 
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FIG. 16. Intrinsic magnetic susceptibility \ °f sample 1 
versus temperature T and fits by the Coqblin-Schrieffer model 
prediction for spins J = 2 and 5/2. The inset shows an 
expanded plot of the data and fits below 40 K. 



V. SUMMARY AND DISCUSSION 

In this paper we have described the synthesis and char- 
acterization of nine LiV204 samples. Our magnetically 
purest samples 1 and 6 clearly showed a broad shal- 
low maximum in the observed magnetic susceptibility 
X obs (T) at T » 16K, with small Curie-like upturns below 
~ 5K. Field-cooled and zero- field-cooled magnetization 
measurements with 77 = 10-100 G did not reveal any ev- 
idence for static spin-glass ordering from 1.8-2 to 50 K 
in any of the seven samples measured. At T > 50 K, 
X° bs (7 1 ) showed local magnetic moment behavior for all 
samples. In sample 2 which showed a larger Curie-like 
upturn in x° b8 (^) a ^ l° w T than in samples 1 and 6, we 
found that liquid-nitrogen quenching reduced the Curie- 
like upturn to a large extent, revealing the broad peak 
in x° bs (7 1 ). However, ice- water quenching and slow-oven 
cooling enhanced the upturn, and the above successful re- 
duction of the upturn by liquid-nitrogen quenching could 
not be reproduced. We analyzed low-T isothermal mag- 
netization versus applied magnetic field M obs (77) data, 
and determined the parameters of the paramagnetic im- 
purities giving rise to the Curie-like upturn in x° bs (7 1 ), 
assuming that a single type of impurity is present. Using 
these parameters, the intrinsic susceptibility x{T) was 
obtained and found to be essentially the same in all sam- 
ples but one (4A) . Surprisingly, the spin Si mv of the para- 
magnetic impurities was found to be large, Si mv = 3/2 
to 4 depending on the sample, suggesting the presence of 
variable amounts of ferromagnetically coupled vanadium 
spin defect clusters of variable size in the samples. 

We tested the localized magnetic moment picture for 
X(T) at T > 50 K using the HTSE prediction for the spin 
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susceptibility of the 5 = 1/2 vanadium sublattice of the 
spinel structure, which yielded C and 9 values similar to 
those reported in the past for LiV204. Using the values of 
the Van Vleck susceptibility obtained from K-\ analyses, 
the Pauli susceptibility contribution to the temperature- 
independent susceptibility xo was derived and found to 
be small, comparable to that of LiTi2 04. The Van Vleck 
formulas for the paramagnetic susceptibility of isolated 
V +3 or V +4 ions or an equal mixture, assuming that 
each V ion is in a cubic CEF, failed to describe the T 
dependence of the observed effective magnetic moment. 
For the high-T "localized moment" region, the observed 
effective moment is in agreement with the spin-only value 
expected for g ~ 2. 

Our attempts to describe the low-T susceptibility data 
in terms of the single-ion Kondo (S = 1/2) and Coqblin- 
Schrieffer {J or S > 1/2) models for isolated mag- 
netic impurities in metals were unsuccessful. These mod- 
els predict that the electronic specific heat coefficient 
7 (T) and the susceptibility x(T) both show maxima for 
J > 3/2.B LiV 2 4 clearly shows a peak in X(T\,ft, 
T w 16 K, but there is no peak in 7 (T) down to 1.2 KM5 
Thus, these theories cannot self-consistently explain the 
results of both measurements, suggesting that there is 
some other mechanism responsible for the heavy-fermion 
behavior and/or that the single-ion picture is inappro- 
priate. It is however intriguing that the experimental 
Wilson ratio Rw ~ 1.7 at IK (Ref. ^) is close to that 
(i?w = 2) predicted for the 5 = 1/2 Kondo model. 

In conventional /-electron heavy fermion compounds, 
local /-electron orbitals and conduction electron states in 
non-/ bands hybridize only weakly, resulting in a many- 
body scattering resonance of the quasiparticles near the 
Fermi energy E-p, a large density of quasiparticle states 
T>(Ep), and hence a large quasiparticle effective mass, 
electronic specific heat coefficient and magnetic spin sus- 
ceptibility at low T. Screening of S — 1/2 local mo- 
ments by conduction electron spins leads to a nonmag- 
netic ground stat e an d a saturating spin susceptibility as 
T — > 0. In Sec. IV, we tested several models for x(T) 
which assume the presence of local magnetic moments in 
L1V2O4 which interact weakly with the conduction elec- 
trons. However, in these models as applied to L1V2O4, 
the itinerant and "localized" electrons must both occupy 
ti g orbitals (or bands derived from these orbitals), rather 
than orbitals of more distinct character. One can imag- 
ine a scenario in which the HF behaviors of LiV204 at 
low T arise in a way similar to that of the /-electron 
HF compounds, if the following conditions are fulfilled: 
(i) the trigonal component of the CEF causes the A\ g 
orbital singlet to lie below the E g orbital doublet; (ii) 
one of the 1.5 d-electrons/V is localized in the ground 
A\ g orbital due to electron-electron correlations^ (iii) 
the remaining 0.5 <i-electron/V occupies the E g doublet 
and is responsible for the metallic character; and (iv) 
the band(s) formed from the E g orbitals hybridize only 
weakly with the A\ g orbital on each V ion. This sce- 
nario involves a kind of orbital ordering; a more general 



discussion of orbital ordering effects is given below. 

The geometric frustration for antiferromagnetic order- 
ing inherent in the V sublattice of LiV204 may be im- 
portant to the mechanism for the observed HF behav- 
iors of this compound at low T. Such frustration in- 
hibits long-range magnetic ordering and enhances quan- 
tum spin ftuctiiations and (short-range) dynamical spin 
ordering £300 These effects have been verified to oc- 
cur in the C15 fee Laves phase intermetallic compound 
(Yo.g7Sco.o3)Mn2, in which the Y and Sc atoms are non- 
magnetic and the Mn atom substructure is identical with 
that of V in L1V2O4. In (Y .97Sco.o3)Mn 2 , Shiga et at 
discovered quantum magnetic moment fluctuations with 
a large amplitude (/x rms = 1.3 /xa/Mn at 8K) in their po- 
larized neutron scattering studyQ They also observed a 
thermally- induced contribution, with fi rms = 1.6/ie/Mn 
at 330 K. Further, Ballou et alSB inferred from their 
inelastic neutron scattering experiments the presence of 
"short-lived 4-site collective spin singlets," thereby sug- 
gesting the possibility of a quantum spin-liquid ground, 
state. A recent theoretical study by Canals and LacroixL3 
by perturbative expansions and exact diagonalization 
of small clusters_of a S = 1/2 (frustrated) pyrochlore 
antiferromagnetLa found a spin-liquid ground state and 
an AF spin correlation length of less than one interatomic 
distance at T = 0. Hence, it is of great interest to carry 
out neutron scattering measurements on L1V2O4 to test 
for similarities and differences in the spin excitation prop- 
erties to those of (Yo.97Sco.o3)Mii2. 

(Yo.97Sco.o3)Mn2 has some similarities in properties 
to those of L1V2O4. No magnetic long-range order^ 
ing was observed above 1.4 K (Refs. p9| , f77|) and 0.02 K,lI 
respectively. Similar to L1V2O4, (Yo.97Sco.o3)Mn2 
shows a large electronic, specific heat coefficient 7(0) « 
160-200 mJ/molK 2 .|2E3 However, the T dependences 
of the susceptibility!— I and 7 (Ref. |7^) are very differ- 
ent from those seen in LiV2 04 and in the heaviest /- 
electron heavy fermion compounds. x obs C^) does not 



show a Curie- Weiss-like be) 
increases with increasing T. 
of T up to at least 6.5 K. 
of Mn with Al, Shiga et 
in (Yo.9sSco.o5)(Mni_ a .Al !B )2 



,vior at high T, but rather 
7(T) is nearly independent 
Replacing a small amount 
found spin-glass ordering 
with x > 0.05.EI The sus- 



ceptibility for x — 0.15 shows a Curie- Weiss-like behav- 
ior above ~ 50 K. The partial removal of the geomet- 
ric frustration upon substitution of Al for Mn might be 
anologous to that in our sample 3 in which structural 
defects evidently ameliorate the frustrated V-V interac- 
tions, leading to spin-glass ordering below ~ 0.8K.tl 

The magnetic properties of materials can be greatly 
influenced when the ground state has orbital degeneracy 
in a high-symmetry structure. Such degenerate ground 
state orbitals can become energetically unstable upon 
cooling. The crystal structure is then deformed to a 
lower symmetry to achieve a lower-energy, non-otbitally- 
degenerate ground state (Jahn- Teller theorem) E3 This 
kind of static orbital ordering accompanied by a struc- 
tural distortion is called the cooperative Jahn- Teller 
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effect The driving force for this effect is the compe- 
tition between the CEF and the lattice energies. Orbital 
ordering may also be caused by spin exchange interac- 
tions in a ma gnet ic system with an orbitally-degenerate 
ground state.c3'E2l The orbital (and charge) degrees of 
freedom may couple with those of the spins in such 
a way that certain occupied orbitals become energeti- 
cally favorable, and consequently the degeneracy is lifted. 
As a result, the exchange interaction-becomes spatially 
anisotropic. For example, Pen et al£3 showed that the 
degenerate ground states in the geometrically frustrated 
V triangular lattice Heisenberg antiferromagnet L1VO2 
can be lifted by a certain static orbital ordering. X- 
ray and neutron diffraction measurements detectecLna 
structural distortions or phase transitions in L1V2O4.&E2I 
However, the presence of orbital degeneracy or near- 
degeneracy suggests that dynamical orbital-charge-spin 
correlations may be important to the physical properties 
of L1V2O4. It is not yet known theoretically whether such 
dynamical correlations can lead to a HF ground state and 
this scenario deserves further study. 

Thus far we and collaborators have experimentally 
demonstrated heavy fermion behaviors of L1V2O4 char- 
acteristic of the heaviest-mass J-eksetron HF systems 
from magiietiKation,H specific heat Jl p. nuclear magnetic 
resonance IjlH thermal expansion,E3o and muon spin 
relaxations measurements. Our magnetization study re- 
ported in this paper was done with high-purity polycrys- 
talline samples from which we have determined the low 
temperature intrinsic susceptibility. Nevertheless, high- 
quality single crystals are desirable to further clarify the 
physical properties. In particular, it is crucial to mea- 
sure the low-T resistivity, the carrier concentration and 
the Fermi surface. In addition, when large crystals be- 
come available, inelastic neutron scattering experiments 
on them will be vital for a deeper understanding of this 
d-electron heavy fermion compound. 

On the theoretical side, new physics may be neces- 
sary to explain the heavy fermion behaviors we observe 
in L1V2O4. We speculate that the geometric frustration 
for antiferromagnetic ordering and/or coupled dynamical 
orbital-charge-spin correlations may contribute to a new 
mechanism leading to a heavy fermion ground state. A 
successful theoretical framework must in any case self- 
consistently explain the radically different properties of 
L1V2O4 and the isostructural superconductor LiTi204. 
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